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TECHNICAL FIELD 

The present invention relates to both orthogonal frequency 
division multiple access OFDMA, to orthogonal Wavelet division 
multiple access OWDMA, to code division multiple access CDMA, and 

10 to multi-scale code division multiple access MS-CDMA, for 
cellular telephone and wireless data communications with data 
rates up to multiple Tl (1.544 Mbps), El (2.048 Mbps), Sonet, 
Ethernet, and higher (>10 Gbps) , and to optical CDMA and optical 
OWDMA. Applications are to wire, wireless local area, wide area, 

15 mobile, point-to-point, and satellite communication networks. 
More specifically the present invention relates to a new and 
novel means for combining MS-CDMA with OFDMA, to a new and novel 
OWDMA which is an orthogonal multi-resolution complex Wavelet 
multiple access generalization of OFDMA, and to a new and novel 

20 means for combining MS-CDMA with OWDMA. This new architecture 
MS-CDMA OFDMA/ OWDMA is an attractive candidate to replace current 
and future OFDMA applications and CDMA applications. 
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BACKGROUND ART 

Current OFDMA art is represented by the applications to the 
wireless cellular communications standards IEEE 802.11a, IEEE 
35 802. llg, IEEE 802.15.3a, IEEE 802.16. OFDMA uses the Fourier 
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transform basis vectors as the orthogonal channelization vectors 
for communications with each basis vector multiplied by a symbol 
which is encoded with a data or pilot signal word. 

The discrete Fourier transform DFT implemented as the fast 
5 Fourier transform FFT is defined in equations (1) . Step 1 defines 
the digital sampling interval T over time, the sampling instants 
t=iT where i is the time index and where the sampling time 1/T is 
at least equal to the complex Nyquist sampling rate to prevent 
spectral foldover. Step 2 is the FFT of the time function Y(i) 

10 and step 3 is the inverse FFT" 1 of the frequency spectrum X(k) 
which derives the original time function Y(i). Step 4 defines the 
NxN orthogonal complex DFT matrix E row vectors E(k) which are 
the DFT harmonic vectors or basis vectors or code vectors or 
channelization vectors. Step 5 is the complex baseband 

15 transmitted signal z(i) for one data block and is equal to the 
inverse FFT transform FFT" 1 of the user symbols x(k). 



Unweighted DFT encoding 
1 Sampling interval of DFT 



(1) 
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t -> 0 



T 



(N-2)T 



(N-l)T 



i -» 0 



N-l 



N-l 



where 1/T > complex Nyquist sample rate 



2 FFT of Y(i) 
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X(k) = FFT [ Y(i) ] 



= ^ E(k, i) Y(i) 



- Si exp(-j2rc ki/N) Y(i) 
3 FFT 1 of X(k) 

Y(i) = FFT _1 [ X(k) ] 
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= N' 1 Z k E* (k,i) X(k) 

= N" 1 Z k exp(j27C ki/N) X(k) 



4 DFT orthogonal harmonic code matrix E 

E = NxN DFT orthogonal DFT matrix 

= [E(k,i)] matrix of elements E(k,i) 
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E(k,i)= exp ( - j 27i ki/N ) harmonic k, time index i 
E(k) = harmonic k basis (code) vector 

= [E(k,0), E(k,l), . . . ,E(k,N-l)] 
EE* = N I 

5 <E(k) ,E* (k' )> = S(k-k') N I 

where I = NxN identify matrix 

E' = complex conjugate transpose of E 
8(k-k')= Dirac delta function 
= 1 for k=k' 
10 =0 otherwise 

5 Transmitted DFT complex baseband signal z(i) for 
one data block 

d(k) = data modulation for user k 

= encoded amplitude A(k) and phase (p(k) 
15 x(k) = transmitted symbol encoded with d(k) 

= A(k) exp(jcp(k) ) 
z(i) = FFT _1 [ x(k) ] 

= N" 1 Z k x(k) E*(k,i) 

20 

OFDMA for IEEE 802. llg in reference [1] is illustrated in 
FIG. 1. The channelization filter h(f) 1 covers a 20 MHz 
frequency band 2 assigned to OFDMA. Plotted is the power spectral 
density PSD=|h(f)| 2 of this channelization filter h(f). A N=64 

25 point fast Fourier transform FFT covers this band 2. Consistant 
with the IEEE specif ication, FIG. 1 refers to the DFT which is 
identical to the analog Fourier transform FT since it is the 
sampled data version of the FT. It is convenient to consider the 
DFT in this invention disclosure as the digital format for the 

30 FFT. The DFT frequency spectrum 3 consists of N=64 equally spaced 
filters 4 across this 20 MHz band. Filter spacing is equal to 
the DFT output rate 1/NT = 0.3125 MHz = 20MHz/64. The DFT time 
pulse p(t) 5 is NT=3.2jxs in length and the total DFT period 6 is 
4.0|xs which allows a 0.8^is guard time for p(t). 
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Throughout this invention disclosure it will be understood 
that the FFT fast algorithm will always be used to implement the 
DFT and the inverse FFT" 1 fast algorithm will always be used to 
implement the inverse DFT" 1 . 
5 OFDMA transmitter encoding of the OFDMA waveform in FIG.l 

is defined in equations (2) . Step 1 lists the parameters and 
definitions and step 2 defines the time domain weighting. Step 3 
is the complex baseband transmitted signal z(i). 



10 OFDMA encoding for transmitter (2) 

1 Parameters and definitions from 1 in equation (1) 
and 

h(i) = 20 MHz band filter impulse response 
p(i) - impulse response of the DFT waveform 
15 = real weighting function in 6 in FIG. 1 

N =64 point DFT 

1/T = 30 MHz sample rate for DFT 

> complex Nyquist rate 
NT = 3.2^s DFT length 
20 1/NT = 0.3125 MHz DFT output rate 

= DFT channel separation 
= DFT tone spacings 
52 channels are used: 4 pilot, 48 data 
12 guard band channels for rolloff of the h(k) 
25 2 Pulse p and band filter h weighting for DFT basis vectors 

p = p©h, convolution of p and h 

= filter transfer function in time domain 

for the combined h and p filters 
- [ p(0), p(l), . . ., p(N-l) ] 
30 3 Transmitted OFDMA. encoded baseband signal z(i) 

for one data block 

d(k) = data modulation for user k 

= encoded amplitude A(k) and phase cp(k) 
x(k) = transmitted symbol encoded with d(k) 
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z(i) 



= A(k) exp(jq>(k) ) 
= FFT _1 [ p(i) x(k) ] 
= N" 1 I k p(i) x(k) E*(k,i) 



5 OFDMA for IEEE 802. llg has the strict orthogonality of the 

DFT(FFT) replaced by cross-correlations between the 48 channel 
tones and other impacts due to the band channelization and pulse 

weighting p©h plus the time errors At and frequency errors Af 
from synchronization errors, multi-path, propagation, and 

10 terminal stresses. These impacts on orthogonality are. low enough 
to allow OFDMA support higher values for the symbol signal-to- 
noise ratio S/N in the detection band that are required for 
higher order symbol modulations. The highest order symbol 
modulation currently is 64 state quadrature amplitude modulation 

15 64QAM corresponding to 6 bits per symbol where 6=log 2 (64) and 
log 2 (o) is the logarithm to the base 2. With rate 3/4 
convolutional coding the highest information rate is 4.5 
bits/symbol = 6x3/4. Required S/N at a BER=1.0e-6 is 
approximately S/N-19 dB. 

20 OFDMA for IEEE 802. llg provides 48 channels over a 20 MHz 

frequency band at a symbol rate equal to 0.25 MHz=l/4.0 jxs from 6 
in FIG. 1 and for a maximum information rate equal to 4.5 
bits/symbol this equals a burst rate of 54 MBps = 4.5x48x0.25. 
Some spread spectrum properties are realized by hopping the 20 

25 MHz band, shuffling the channel assignments over the 48 available 
channels for a user assigned to several channels in order to 
spread his transmissions over the band, and for "flash" ODFMA by 
a random hopping of each user channel across the 48 available 
channels within the band. 

30 OFDMA receiver decoding of the OFDMA waveform in FIG.l is 

defined in equations (3) for the receiver. Step 1 identifies the 

parameters and definitions and defines the convolution R(i)= p©p 
of the time domain weighting p with itself in the matched filter 
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receiver. Step 2 derives the estimates of the transmitted symbols 
by implementing matched filter detection in the receiver. 

OFDMA decoding for receiver (3) 
5 1 Parameters and definitions are defined in 1,2 in 

equations (2) together with 

R(i) = convolution of p with p 

= p©p 

2 OFDMA decoding for one data block derives estimates 
10 x(k) of x(k) from the receiver estimates z(i) of z(i) 

x(k) = FFT [ z(i)©p ] 

= Zi z(i)©^E(k,i) 

15 Current CDMA spread spectrum art is illustrated by the 

waveform in FIG. 2 which describes the waveform for full band 
CDMA communications over the band B 9 which is the output 
range 8 of the band filter h(f) 7. The CDMA chip rate 1/T C 12 
is less than the available frequency band B to allow the chip 

20 frequency spectrum p(f) 10,11 to roll off. As defined 9 in 
FIG. 2 the band B and chip rate are related by equation B=(l + 
a) /T c where a is the bandwidth expansion factor close to a=0.25 
for high performance communications. Frequency spectrum p(f) 10 
for the CDMA communications is essentially equal to the Fourier 

25 transform of the CDMA time pulse p(t) 13. The representative 
time pulse p(t) is a square-root raised cosine pulse which is 
typically used for high performance communications to obtain a 
reasonably flat spectrum with a sharp rolloff at the edges to 
enable the chip rate 1/T C to be reasonably close to the available 

30 frequency band B. 

Chip rate 1/T C is the CDMA total symbol rate. The users 
could be at different data rates but this and other architectural 
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variations do not limit the scope of this invention. Power is 
uniformly spread over the CDMA pulse waveform spectrum p(f) . 

It is self evident to anyone skilled in the CDMA 
communications art that these communications mode assumptions 
5 are both reasonable and representative of the current CDMA art 
and do not limit the applicability of this invention. 

CDMA encoding of the waveform in FIG. 2 for the transmitter 
is defined in equations (4) . Steps 1,2 define the CDMA 
transmission and parameters. Step 3 defines the user symbol 
10 x(u). Step 4 is the set of Walsh orthogonal channelization codes 
w(u) and step 5 is the pseudo-random PN covering or spreading 
code. Step 6 defines the complex baseband signal z(t) as the 
waveform p(t-nT c ) multiplied by the encoded sum over u and n. 



15 CDMA encoding for transmitter (4) 

1 CDMA transmission 

chip 0 chip 1 chip N c -1 

t t t 

t -» 0 T c (N c -l)Tc 

20 n -> 0 1 N c -1 

p ^ p(t) p(t-T c ) p(t-(N c -l)T c ) 



where p(t-nT c ) = p(t-nT c ) ©h(t) 
25 = convolution of p(t-nT c ) and h 

= filter transfer function in time domain 
for the combined p(t-nT c ) and h filters 

2 Parameters and definitions 

N c = Number of users and orthogonal code chips 
30 T c = CDMA chip length or repetition interval 

1/N C T C = User symbol rate 

3 User complex signal x(i) 

d(u) = data modulation for user u 

= encoded amplitude A(u) and phase (p(u) 
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x(u) = transmitted symbol encoded with d(u) 
= A(u) exp(jq>(u) ) 

4 Walsh orthogonal channel izaat ion code matrix W 

W = Code matrix, N c rows of N c code vectors 

= [ W(k,n) ] matrix of elements C(k,n) 

W(u,n)= +/-!/ chip n of vector u 

W(u) = code vector u, row k of W 

= [W(u,0), W(u,l), . . . ,W(u,N c -l)] 

5 PN covering (spreading) code P(n) for chip n 

P (n) = exp ( jq> (n) ) 

6 Transmitted CDMA complex baseband signal z(t) 

z(t) = Nc" 1 Z u S n P(n)W(u,n)x(u) p(t-nT c ) 

CDMA decoding of the waveform in FIG. 2 for the transmitter 
is defined in equations (5) . Step 1 defines the convolution 
R(n,n-n') of the CDMA pulse waveform with itself in the matched 
filter receiver. Steps 2,3 are the Walsh, PN decoding 
properties. Step 4 uses the matched filter detection theorem to 
derive the estimates of the transmitted symbols. 

CDMA decoding for receiver (5) 

1 Parameters and definitions are defined in 

1,2,3 in equations (4) together with 

R(n,n-n') = convolution of p(t-nT c ) with 

p(t-n'T c ) evaluated at the receiver 
detection times t=nT c 
= R(n) 8(n-n' ) 
= R(n) for n=n' 
= 0 otherewise 

2 Walsh decoding of channelizaation codes W(k) 

WW* = N c I 

where I = N c x N c identify matrix 

W* = conjugate transpose of W 
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<W(k) ,W(k' ) > = N 5(k-k' ) 
3 PN decoding 

P(n) P* (n) = 1 for all n 
where P* = complex conjugate of P 
5 4 CDMA, decoding 

x(k) = Z n P*(n)W*(k,n) z(t)© £ (t-nT c ) 

It should be obvious to anyone skilled in the 
communications art that these example implementation algorithms 

10 in equations (1) , (2) , (3) , (4) , (5) clearly define the fundamental 
OFDMA and CDMA signal processing relevant to this invention 
disclosure and it is obvious that this example is representative 
of the other possible signal processing approaches. 

For cellular applications the encoding algorithms for the 

15 transmitter describe the implementation of OFDMA and of CDMA 
encoding and are the transmission signal processing applicable to 
this invention for both the hub and user terminals, and the 
decoding algorithms for the receiver describes the corresponding 
OFDMA and CDMA receiving signal processing for the hub and user 

20 terminals for applicability to this invention. 

For optical communications applications the microwave 
processing at the front end of both the transmitter and the 
receiver is replaced by the optical processing which performs the 
complex modulation for the optical laser transmission in the 

25 transmitter and which performs the optical laser receiving 
function of the microwave processing to recover the complex 
baseband received signal with the remainder of the signal 
processing functionally the same for the OFDMA and for the CDMA 
encoding transmitter and functionally the same as described for 

30 the OFDMA and CDMA receiving signal processing receiver. 
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SUMMARY OF INVENTION 

This invention introduces the new OWDMA communications 
technology which implements orthogonal multi-resolution complex 
Wavelet division multiple access and is a multi-resolution 
5 complex Wavelet generalization of OFDMA, introduces the new 
application of the multi-scale code division multiple access MS- 
CDMA architecture which integrates MS-CDMA with OFDMA and with 
OWDMA, and introduces the variable gain control over frequency. 
The new OWDMA forms a uniform set of contiguous orthogonal 

10 filters across the available frequency band in one of the 
numerous available architectural options and which implements a 
polyphase filter bank across the available frequency band with 
the basic property that the filters are orthogonal. Each filter 
defines a OWDMA channel for communications. Similar to OFDMA, 

15 the symbol rate within each channel is equal to the channel 
separation 1/NT where N is the number of channels, 1/T is the 
frequency band, and NT is the symbol-to-symbol separation and 
which is made possible by using multi-resolution complex Wavelet 
channelization filters developed in reference [2] . OWDMA 

20 filters are orthogonal in frequency which means their frequency 
spectrums are non-overlapping, and they have flat spectrums 
across the channel. 

OWDMA orthogonality only requires frequency synchronization 
whereas OFDMA orthogonality requires both time and frequency 

25 synchronization which means that time synchronization errors on 
the user communication channels for the return links to the hub 
or access point for cellular communications does not degrade the 
orthogonal separation between the channels as they will for 
OFDMA. 

30 Sensitivity to frequency synchronization errors is less for 

OWDMA primarily because of the non-overlapping of the frequency 
spectrums. The overall tolerance to user-to-user imbalances on 
the return communications channels is better for OWDMA primarily 
because their design keeps the frequency spectrums from 

35 overlapping in the presence of real operational conditions with 



synchronization errors. This allows the forward communications 
link to support a power imbalances between channels to mitigate 
differing ranges and path losses at the user antennas. 

Application of multi-scale code division multiple access 
5 MS-CDMA in reference [2] provides a means to implement the new 
hybrid-Walsh orthogonal CDMA codes in reference [4] over the 
OFDMA/OWDMA channels by spreading the CDMA within each channel 
and over all of the channels such that each user can be spread 
over the complete band. This keeps the chip rate equal to 1/NT 

10 while maintaining the spreading over the fullband 1/T and allows 
the band Tx power to be independently controlled in frequency. 
These two CDMA scales 1/NT, 1/T are generated by MS-CDMA in 
combination with the OFDMA/OWDMA channelization filter banks. The 
1/T scale is to combat fading and interference similar to the 

15 current CDMA, and the 1/NT scale is for acquisition, 
synchronization, and protection against multi-path and provides 
the flexibility for band power control to provide a frequency 
diversity communications improvement . 

Variable control over the frequency B can be implemented by 

20 partitioning the CDMA over the channels into separate groups and 
assigning an independent power level to each group of channels. 
OWDMA readily supports differences in power levels between 
adjacent channels. Power control is desirable to support 
differences in quality of service SoC, range losses , and path 

25 loss. 

It should be obvious to one familiar with the CDMA 
communications art that the number of scales could be a larger 
number than the two used in this invention disclosure and the 
multi-resolution complex Wavelet design for OWDMA supports the 

30 partitioning of the frequency band 1/T into several frequency 
scales for the channelization filters, supports individual multi- 
resolution complex Wavelet packets for communications in the 
time-frequency domain and which can be integrated into the MS- 
CDMA, and supports separate and segmented communications bands 

35 simultaneously . 
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BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a description of the OFDMA waveform in the 
5 frequency and time domains over a 20 MHz band for the IEEE 
802. llg standard. 

FIG. 2 is a description of the CDMA waveform in the 
frequency and time domains over a band B. 

FIG. 3 is a description of the OWDMA waveform in the 
10 frequency and time domains over a band B. 

FIG. 4 is a description of the design requirements on the 
multi-resolution complex Wavelet PSD= 1 1|/ (f ) I 2 . 

FIG. 5 is a representative baseband frequency response for 
the multi-resolution complex Wavelet filter and the square-root 
15 cosine filters for the bandwidth expansion factors =0.22, 0.40. 

FIG. 6 is a representative implementation of MS-CDMA 
OFDMA/ OWDMA encoding for transmitter. 

FIG. 7 is a representative implementation block diagram for 
the MS-CDMA OFDMA/ OWDMA transmitter. 
20 FIG. 8 is a representative implementation of MS-CDMA 

OFDMA/OWDMA encoding for the receiver. 

FIG. 9 is a representative implementation block diagram for 
the MS-CDMA OFDMA/OWDMA receiver. 



30 
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DISCLOSURE OF INVENTION 

This invention introduces the new orthogonal multi- 
resolution complex Wavelet division multiple access OWDMA 
communications technology which is a multi-resolution complex 
5 Wavelet generalization of OFDMA, introduces the new multi-scale 
code division multiple access MS-CDMA architecture to integrate 
MS-CDMA with OFDMA, and introduces the new MS-CDMA architecture 
to integrate MS-CDMA with OFDMA, 

OWDMA uses the multi-resolution complex Wavelet waveform 

10 developed in reference [2] to generate multi-rate orthogonal 
filter banks and waveforms to support communications at the 
critical symbol rate equal to a combined 1/T symbols per second 
for a 1/T Hz frequency band. The critical symbol rate is the 
Nyquist sample rate and corresponds to no excess bandwidth oc=0 

15 in the equation B=(l+a)/T which can be rewritten as the 
bandwidth-time product BT=(l+a) that relates filter bandwidth B 
and the symbol rate 1/T for communications supported by this 
filter. For the OFDMA filter in FIG. 1 it is observed 
BT=l+a=64/52=l . 33 corresponding to a=0.33 since 52 of the 64 

20 DFT filters generated are used and for the CDMA filter in FIG. 2 
it is observed that the bandwidth-time product is about BT=1.25 
corresponding to a=0.25 for the 3G CDMA cellular communications 
systems . 

Multi-resolution complex Wavelet design algorithms were 
25 developed in reference [2] as a means to design polyphase 
multirate filters, quadrature mirror filters (QMF) , perfect 
reconstruction filters, Wavelet iterated filter banks, and 
Wavelet tiling of the time-frequency t-f space. Prior to the 
invention of multi-resolution complex Wavelet design algorithms, 
30 theorectical studies had not yielded useful realizable filters 
for system applications implementing these architectures as 
summarized by the digital filtering and polyphase research in 
references [5]. [6] and the Wavelet research as summarized in the 
references [7] , [8] , [9] ,10] , [pi] . 
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Multi-resolution complex Wavelets defined in equations (6), 
(7) , (8) expand the Wavelet analytical formulation to include a 
frequency variable which specifies the center frequency of our 
new Wavelets and performs a frequency translation of the real 
5 mother Wavelet to this center frequency. Currrently, Wavelets are 
real functions of the scale and translation parameters. Multi- 
resolution complex Wavelets are functions of these parameters 
plus the frequency variable. This new concept of frequency as an 
additional parameter provides an added degree of flexibility and 

10 together with the Fourier domain design approach provide an 
entirely new means for deriving these new waveforms as 
generalization of the traditional real Wavelets and their 
modification to include the added frequency parameter. With 
frequency translation as the additional parameters the analytical 

15 formulation of these new waveforms as a function of the baseband 
or mother waveform centered at dc ( dc refers to the origin f=0 
of the frequency space) is given in equations (6) . 

Equations (6) introduce the new multi-resolution complex 
Wavelets. Step 1 in equations (6) is the definition for the 

20 continuous real Wavelet over the time-frequency t-f space from 
reference [7], [9], [12] where i|/(t) is the Wavelet which is a 
waveform of finite extent in time t and frequency f over the t-f 
space, Wavelet parameters a,b are the Wavelet dilation and 
translation respectively or equivalently are the scale and shift. 

25 The iy(t) without the indices a,b is the mother Wavelet which is a 
real and symmetric localized function in the t-f space used to 
generate the doubly indexed Wavelet iy(t|a,b) where \|/(t|a,b) reads 
the V)/ is a function of time t for the parameters a,b. Scale 
i r 1/2 

factor a has been chosen to keep the norm of the Wavelet 
30 invariant under the parameter change a,b. Norm is the square root 
of the energy of the Wavelet response. The Wavelets \|/(t|a,b) and 
v|/(t) are localized functions in the t-f space which means that 
both their time and frequency lengths are bounded. 
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Step 2 in equations (6) defines the digital Wavelet which 
is a Wavelet over the digital t-f space corresponding to 
digitization of the t-f space at the 1/T analog-to-digital A/D 
rate where T is the interval between digital samples and "i' is 
5 the index over the digital samples. Time x> t" for the continuous 
Wavelet in 1 is replaced by the equivalent digital sample number 
"i" corresponding to t=iT at sample i. Wavelets in digital t- 
f space have an orthogonal basis that is obtained by restricting 
the choice of the parameters a,b to the values a=2" p , b=qN2 p 
10 where p,q are the new scale and translation parameters and N is 



New multi -resolution complex Wavelets 



(6) 



1 Current continuous Wavelet as a function of the 



mother Wavelet at dc 
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y(t|a,b) = Ia|" 1/2 (t-b) /a) 
where a,b are the dilation, shift parameters 



2 Current digital real Wavelet 
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Digital Wavelet shift q and scale p 
a = 2 P 
b = q N 2 P 



Digital Wavelet as a function of mother Wavelet 
y (i|p,q) = 2~ p/2 i)/(2- p i-qN) 



where N = number of samples over the Wavelet 
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spacing or repetition interval 
T s = NT wavelet spacing 



i = digital sampling index 

T = digital sampling interval 

t = iT time index at digital sample i 



3 Our new digital multi-resoltion complex Wavelet 
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is a function of the mother Wavelet at dc 



3. a) Generalized complex format 

i|/(i|p,q,k) = 2" p/2 vj/(2- p i-qN) exp [ j 27if c (p, k) 2" p iT] 



3.b) Multi-rate filter complex format 



i|/(i|p,q,k) = 2" p/2 M/(2" p i-qN)E* (k2 p , 2" p i) 



4 Orthogonality equations for our new digital 
complex Wavelet 

EiV(i|p,q f k) if/*(ilp',q',k') =N if f p=p' , q=q' , k=k' 

= 0 otherwise 

5 New digital complex Wavelet for a uniform filter bank 

as a function of the mother Wavelet at dc 
v|/(i|p=0,q,k) = \|/(i-qN) E*(k,i) 

the spacing or repetition interval T S =NT of the Wavelets (which 
from a communications viewpoint are symbols) at the same scale p. 
Wavelets at p,q are related to the mother Wavelet by the equation 
in 2 where the mother Wavelet is a real and even function of the 
sample coordinates and which follows directly from 1 for the 
continuous t-f space. 

Steps 3,4 in equations (6) define our new Wavelet in the 
digital t-f space and their orthogonality properties. Our new 
Wavelets are complex generalizations of Wavelets in t-f space 
which enable them to be useful for communications and radar 
applications. This generalization is accomplished 1) by the 
addition of a frequency translation parameter k which controls 
the frequency offset of the Wavelet, 2) by generalizing the 
Wavelet weighted orthonofmality properties in step 4 to apply to 
waveforms over the time translation q and over the scales p with 
the inclusion of the frequency translation and where is the 
complex conjugate of i|/ and is required in the orthogonality 
equations since the multi-resolution complex Wavelet vj; becomes 
complex with the addition of harmonic k, and 3) by their 
characterization and design in the Fourier domain. The frequency 
parameter k controls the frequency translation exp [ j27tf c (p, k) 2" 
p iT] in 3 for the generalized format and the frequency 
translation E*(k2 p ,2" p i) in 3 for the multi-rate filter format. 
With this frequency translation the analytical formulation 3 of 
these new Wavelets is given as a function of the baseband or 
mother waveform centered at dc corresponding to k=0. Purpose of 
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the frequency index k is to identify the center frequencies of 
the waveforms at the scale p and time translation q in the t-f 
space. The generalized center frequency f c (p,k) of the frequency 
translated dc waveform at scale p and frequency index reduces to 
5 f c (P/k)= k2 p /NT for application to multi-rate filters. 

Step 5 in equations (6) defines the equations for our new 
Wavelet for a uniform polyphase filter bank which is one of 
several OWDMA candidate architectures. Our new Wavelets in 5 are 
the impulse responses of the corresponding digital symbols 

10 encoded with digital data for transmission which is the synthesis 
filter bank in polyphase theory, and are the filter bank 
detection impulse responses for recovery of the transmitted 
digital symbols by the analysis filter bank in polyphase theory 
assuming matched filter detection. The digital filters are 

15 observed to be the DFT's of the mother wavelets similar to the 
construction of the OWDMA in equations 3 in equations (2) upon 
replacing the OWDMA pulse waveform p with the mother Wavelet 
H/(i-qN) in 5 and setting q=0 corresponding to the symbol at t=0. 
Equations (7) derive the multi-resolution complex Wavelet 

20 as a function of the Fourier design coordinates. Design 
algorithms provide a means to design the mother Wavelet in the 
Fourier frequency domain to fit the communications and radar 
specifications. From this mother Wavelet, the Wavelets at the 
appropriate scales {p,q,k} are easily found as demonstrated in 

25 3,5 in equations (6) . Design in the frequency domain means the 
design coordinates specifying the Wavelet are Fourier frequency 
harmonics or coordinates. Step 1 list the parameters and 
coordinates and step 2 define the Fourier harmonic frequency 
design coordinate. 

30 Step 3 is the DFT representation from equations 3 in 

equations (1) of the real mother Wavelet v|/(i) as a function of 
the DFT Fourier harmonic coefficients {\|/(k' )}. Step 4 derives 
the equation for the multi-resolution complex Wavelet as a 
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function of the DFT Fourier harmonic coefficients by substituting 
3 in equations (7) into 3 in equations (6) . 

Digital multi-resolution complex Wavelet is a (7) 
function of the Fourier harmonic design coordinates 

1 Parameters and coordinates 

\\f = v|/(i) mother Wavelet 

= \|/(Up=0,q=0,k=0) 
N' = length of the mother Wavelet \|/(i) 

= NL+1 

N = sampling interval of \j/ 

= spacing of \\j for orthogonality 
L = length of V|/ in units of M 

T = spacing of digital samples 

T s = N'T length of \\f in seconds 

2 Fourier harmonic frequency design coordinates 

i|/(k') = value of the Fourier design harmonic k'/N'T 

= Fourier harmonic design coordinate value 
{k'} = set of harmonic design coordinates 
for k' = 0, . . . 

corresponding to harmonics {k'/N'T} 

3 Multi-resolution real mother Wavelet definition 

V|/(i) = Z k . V|/(k') E*(k',i) 

4 Multi-resolution complex Wavelet definition 

i|/(i|p,q,k) 

= 2" p/2 £ k ,i|/(k') E* (k' /L+k2 p , 2~ p (i-qN) ) 

Several fundamental properties follow directly from the 
frequency design approach in equations 4 in equations (7) . It is 
demonstrated in reference [2] that our multi-resolution complex 
Wavelets are implemented with our design in the Fourier domain 
and our multi-resolution complex Wavelet design remains invariant 
under scale changes. It is demonstrated in reference [2] that the 
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Fourier frequency domain design in equations 3 in equations (7) 
remains invariant for all parameter changes and in particular for 
all scale changes. 

Multi-resolution complex Wavelet design algorithms are 

5 illustrated by a representative least-squares LS design for a 
OWDMA polyphase uniform filter bank in FIG. 3. There are two 
categories of LS algorithms and these are the eigenvalue and the 
gradient search that respectively can be reduced to algorithms 
equivalent to the original eigenvalue [13] and Remez-exchange 
10 [14] waveform design algorithms for application to a uniform 
filter bank. 

FIG. 3 is an example of an orthogonal multi-resolution 
complex Wavelet division multiple access OWDMA Wavelet polyphase 
filter bank over the band B by setting N in 5 in equations (6) 

15 equal to the N in 17 in FIG. 3. The band channelization filter 
is a roofing filter h(f) 14 that covers the B frequency band 15 
assigned to OWDMA. Unlike its use in FIG. 1 for the OWDMA with 
B=20 MHz this filter is intended to be flat over the frequency 
band B of interest and with a bandwidth 18 equal to 1/T > B. 

20 Plotted is the power spectral density PSD=|h(f)| 2 of this 
channelization filter h(f). The h(f) output is digitized at the 
sample rate 1/T to form the OWDMA multi-resolution complex 
Wavelet polyphase filter bank 16 which are uniformly spaced 17 
at 1/NT Hz over the 1/T frequency band 1/T. This digital sample 

25 rate 1/T is sufficiently large to allow the use of the shaded 
OWDMA filters in 17 in FIG. 3 for communications over the band 
B 18 with no excess bandwidth a=0 unlike the OFDMA in FIG. 1 
which has a=0.33 and CDMA in FIG. 2 with <x=0.25. 

This OWDMA polyphase filter bank is ideally decimated which 

30 means the filter output multi-resolution complex Wavelet sample 
rate 1/T S is equal to the channel-to-channel spacing T S =NT 
equivalent to stating that there is no excess bandwidth a=0 
within the filter bank. A representative 19 OWDMA multi- 
resolution complex Wavelet \\f for L=8 is plotted in 20 as a 
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function of the time offset expressed in units of the multi- 
resolution complex Wavelet spacing NT. This i|/ was designed by 
the eigenvalue category of LS design algorithms. Our design for 
this topology is immediately applicable to an arbitrary set of 
5 multi-resolution OWDMA filters through the scaling equations 3,4 
in equatioins (7) which gives the design of our Wavelet at 
arbitrary scales in terms of our design of the mother Wavelet. 

Multi-scale mother Wavelet \\f design for the OWDMA polyphase 
filter bank in FIG. 3 starts by using the frequency design 

10 template in FIG. 4 to construct the LS error metrics as a 
function of the DFT frequency design coordinates {\\f (k' ) in 2 in 
equations (7) which define the mother Wavelet in 3 in equations 
(7) . Next minimization search algorithms are developed and used 
to find the values of {i}/(k') that minimize the weighted sum of 

15 these LS design error metrics equal to the cost function J in 
equations (8) for the LS design. Minimizing J with respect to 
the {i]/(k') gives the design values of {i|/(k'} for constructing 
the mother Wavelet in equations 3 in equations (7) . 

Equations (8) define the LS cost function J constructed 

20 with the LS metrics with the aid of the frequency design template 
in FIG. 4. Step 1 define the LS error metrics which are the 
passband error metric fi(l) , stopband error metric \i{2) , 
quadrature mirror filter QMF error metric (J.(3) , intersymbol 
interference error metric jx { 4 ) , and the adjacent channel 

25 interference error metric n(5) . 

LS cost function J for designing V(/(i) (8) 
1 LS metrics 

|x(l) = passband LS error metric measures the LS 
30 error of the passband ripple 24 in FIG. 4 

\i(2) = stopband LS error metric measures the 

stopband atteunation 28 in FIG. 4 
\x (3) = quadrature mirror filter QMF LS error 
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metric measures the flatness of the sum 
of two contiguous filters over the 
deadband 31 in FIG. 4 

|u(4) = orthogonality LS error metric measures the 
5 intersymbol interference ISI between 

overlapping Wavelets of different wymbol 
[i(5) = orthogonality LS error metric measures the 
adjacent channel interference ACI from 
nearest neighbor channels 
10 2 Metric weighting 

w(n) = metric weight for error metrics n=l,2,3,4,5 in 1 
> 0 

such that I n w(n) =1 
3 Cost function J 

15 J = S n w(n) n(n) 

FIG. 4 is the frequency design template for the power 
spectral density PSD= | \\t (f ) | 2 of the multi-resolution complex 
Wavelet and defines the parameters of interest for the passband 

20 metric |i(l) and stopband metric |u,(l) in step 1 in equations (8) . 
Passband 21 of the wavelet PSD=|\|/(f) | 2 is centered at dc (f=0) 
since we are designing the mother Wavelet, and extends over the 
frequency range (Dp extending from -C0p/2 to +G)p/2 22 in units of 
the radian frequency variable co=27ifT 23 where T is the 

25 digital sampling interval defined in FIG. 3. The frequency space 
extends over the range of Af=-1/2T to Af=+1/2T which is the 
frequency range in FIG 3 and the mother Wavelet is at the center 
of the frequency band. Quality of the PSD= | \\i (f ) | 2 over the 
passband is expressed by the passband ripple 24. Stopband 25 

30 starts at the edge 26 of the passbands of the adjacent channels 
+ /-G> a /2 26 and extends to the edge of the frequency band co=+/- 
n 27 respectively. Stopband attenuation 28 at +/-(D a /2 

measures the PSD= | \\f (f ) | 2 isolation between the edge of the 
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passband for the mother Wavelet and the start of the passband for 
the adjacent Wavelet channel centered at +/-cg s 29. Rolloff 30 
of the stopband is required to mitigate the spillover of the 
Wavelet channels other than the adjacent Wavelet channels , onto 
5 the mother Wavelet channel. Deadband or transition band 31 is 
the interval between the passbands of contiguous Wavelet 
channels, and is illustrated in FIG. 4 by the interval from co p /2 
to G> a /2 between the mother Wavelet channel and the adjacent 
Wavelet channel at oo a . Waveform sample rate co s 32 is the 

10 waveform repetition rate. For the LS example algorithm, the 
waveform sample rate is equal to the channel-to-channel spacing 
for zero excess bandwidth. Therefore, 1/T S = co s /27iT=l/MT which 
can be solved to give <fl s =27i/M for the radian frequency sampling 
rate of the filter bank which is identical to the Wavelet 

15 repetition rate. 

Equations (8) step 1 QMF LS error metric \i(3) expresses the 
requirement on the deadband that the PSD' s from the contiguous 
channels in FIG. 4 add to unity across the deadband 31 [G) p , 
co s ] in FIG. 4 in order that the Wavelets be QMF filters. 

20 Equations (8) step 1 ISI and ACI error metrics |i(4),|x(5) are 

orthonormality metrics that measure how close we are able to 
designing the set of Wavelets to be orthonormal over the t-f 
space, with the closeness given by the ISI error metric |a ( 4 ) and 
the ACI error metric n(5) . ISI is the non-orthogonality error 

25 between Wavelets within the same channel separated by multiples 
of the sampling interval 1/MT seconds where T is the sample time 
and M is the interval of contiguous samples. ACI is the non- 
orthogonality error between between Wavelets within a channel and 
the Wavelets in adjacent Wavelet channels at the same sample time 

30 and at sample times separated by multiples of the sample 
interval. As observed as noise contributions within each sample 
in a given channel, the ISI is the noise contribution due to the 
other received Wavelets at the different timing offsets 
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corresponding to multiples of the sampling interval. Likewise, 
the ACI is the noise contribution due to the other Wavelets in 
adjacent Wavelet channels at the same sampling time and at 
multiples of the sampling interval. 
5 Equations (8) step 1 ISI and ACI errors are fundamentally 

caused by different mechanisms and therefore have separate 
metrics and weights to specify their relative importance to the 
overall sum of the LS metrics. ISI is a measure of the non- 
orthogonality between the stream of Wavelets within a channel as 

10 per the construction in FIG. 4. On the other hand, ACI is a 
measure of the non-orthogonality between the Wavelets within a 
channel and the other Wavelets in adjacent channels. This means 
the stopband performance metric has a significant impact on the 
ACI due to the sharp rolloff in frequency of the adjacent 

15 channel, and the ACI metric is then a measure of the residual 
non-orthogonality due to the inability of the stopband rolloff in 
frequency from completely eliminating the ACI errors. 

Equations (8) step 2 defines the weights of the LS error 
metrics when summed to yield the cost function J. These weights 

20 are real and normalized to sum to unity. They have proven to be 
helpful in the Wavelet design to emphasize the relative 
importance of the individual error metric contributions to J. 

Equations (8) step 3 defines the cost function J as the 
weighted sum of the LS error metrics and which is minimized with 

25 respect to the DFT frequency design harmonics {\j/(k')} to select 
the best LS choice for the {\|/(k')> to design the mother Wavelet 
in 3 in equations (7) and the channel Wavelets by frequency 
translation in 5 in equations (6) . 

Multi-scale mother Wavelet frequency response in FIG. 5 is 

30 evaluated by implementing the LS design algorithms in reference 
[2] that minimize the J in 3 in equations (8) to find the best 
set of DFT frequency design coordinates { vj/ ( k' ) } which give the 
mother Wavelet 3 in equations (7) . FIG. 5 plots the PSD 
frequency response for the multi-scale mother Wavelet and the 
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square-root (sq-rt) raised-cosine (r-c) waveforms with excess 
bandwidth a= 0.22, 0.40 which waveformns are extensively used for 
the third generation 3G CDMA communications as well as for other 
communications. Plotted are the measured PSD in dB 33 versus 
5 the frequency offset 34 from dc expressed in units of the symbol 
rate Plotted against the normalized frequency offset are the 
mother PSD for the new waveform 35, the sq-rt r-c with a=0.22 
36, the sq-rt r-c with a=0.40 37, passband template for the dc 
filter and the stopband template for the contiguous filter, and 
10 the Nyquist rate equal to the multi-resolution complex Wavelet 
rate = 0.5 in units of the symbol rate. It is believed that the 
multi-resolution complex Wavelets can be designed as a filter 
with better performance parameters than possible with any other 
known algorithm. 

15 OWDMA encoding for the transmitter is defined in 

equations (9) . Step 1 lists parameters and definitions and step 
2 defines the transmitted OWDMA encoded baseband signal z(i) for 
contiguous data blocks and where the symbol offsets A account for 
symbol overlaps over the symbol q data block interval within each 

20 channel. 

OWDMA encoding for transmitter (9) 
1 Parameters and definitions 

h(i) = roofing filter impulse response in time 
25 for h(f) in 15 in FIG. 3 

v|/(i) = \|/ (i | p=0, q=0, k=0) mother Wavelet in 3 
in equations (7) 
= baseband or dc k=0 Wavelet at p=0,q=0 
N = number of OWDMA filters over the 1/T 

30 frequency band 

1/T = digital sample rate for OWDMA 

> complex Nyquist rate for roofing filter h(f) 
NT = OWDMA Wavelet spacing 
1/NT = OWDMA Wavelet output rate 
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= OWDMA Wavelet channel separation 
NBT = channels used for data and pilot 

where B is the frequency band in FIG. 3 
N(l-BT) = guard band channels for rolloff 
5 of the h(f) 

d(k,q) = data modulation for user k for data block q 

= encoded amplitude A(k|q) and phase (p(k|q) 
x(k|q) = transmitted symbol encoded with d(k) 
= A(k|q) exp(jq>(k|q) ) 
10 Assume the h(f) is flat over the passband for 

both Tx and Rx and can be neglected 
2 Transmitted OWDMA encoded baseband signal z (i) 
Time index field referenced to q=0 
i = 0, +/-2, . . . 

15 = digital sampling time index 

i 0 = index over a i|/(i) spacing interval 

= 0,1,2,. . . . ,N-1 
A = index over the i|/(i) range in units 
of the \\f(i) spacing interval N 
20 = (L/2-1) , . . .,-l,0,+l, . . .,+(L/2-l) 

Complex baseband signal 

z(iolq) = complex baseband signal over i 0 
for data block q 
= FWT[ x(k|q+A) ] 
25 = N' 1 Z A E k x(k|q+A)\]/(i|q+A,k) 

= N" 1 E A Z k x(k|q+A)i|/(i 0 +AN)E* (k,i 0 ) 

z(i) = S q z(i 0 |q) 
where FWT is the fast multi-resolution complex 
Wavelet transform 
30 3 FWT algorithm for OWDMA encoding in the transmitter 

3. a) FWT pre-calculation FFT" 1 

Mi 0 ,q+A) = N" 1 Z k x(k|q+A) E*(k,i 0 ) 

= FFT' 1 [ N" 1 E k x(k|q+A) E*(k,i 0 ) ] 
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3.b) FWT post-sum 

z(iolq) = £ A H/(i 0 +AN) Miolq+A) 



4 Computational complexity of fast algorithm 

5 Real multiply rate R M 

R M T = 2 log 2 (N) + 2 L 
Real add rate R a 

R A T = 3 log 2 (N) + 2 L 

Step 3 is the new fast FWT algorithm in this invention 
disclosure for the transmitted OWDMA which consists of the pre- 
calculation FFT -1 in sub-step 3. a followed by a post-sum in sub- 
step 3.b. 

Step 4 evaluates the real multiply complexity metric R M T 
and real add computational complexity metrics R A T in terms of 
multiplies/adds per digital sample for the fast algorithm in step 
3. The first term in these metrics is the complexity of the FFT" 1 
for a base 2 and the second term is the complexity of extending 
the multi-resolution complex Wavelet waveform over L of the 
symbol intervals. 

OWDMA decoding for the receiver is defined in 
equations (10) . Step 1 refers to 1,2 in equations (9) for the 
parameters and definitions and defines the OWDMA filtering 
Wavelet. Step 2 demonstrates Wavelet orthogonality. Estimates of 
the transmitted symbols in step 3 are equal to the FWT" 1 of the 
received baseband signal. 

OWDMA decoding for receiver (10) 
1 Parameters and definitions are defined in 1,2 in 
30 equations (9) together with 

OWDMA filtering wavelet from 5 in equations (6) 

V*(i|q,k) = n/*(ilp=0,q,k) 

= i|/(i-qN)E(k,i 0 ) 
Assume the h(f) is flat over the passband for 
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both Tx and Rx and can be neglected 
2 Multi-resolution complex Wavelet orthogonality from 4 in 
equations (6) 

k,k' orthogonality 

5 EiV(i|q,k) i|/*(i|q,k') 

= Iii|/ 2 (i-qN)E*(k,i)E(k',i) 
= E i0 Sa [2 Aq v|/ 2 (i 0 +A)N)] E*(k,i 0 )E(k',i 0 ) 
= E i0 [ 1 ]E*(k # i 0 ) E(k' ; i 0 ) 
= N 8(k-k' ) 

10 using the completeness property of the 

multi-resolution complex Wavelet from reference [2] 
[Z A i|/ 2 (i 0 + A) N) ] = 1 for all i 0 
q,k and q' ,k' orthogonality proven in reference 
[2] and restated in 4 in equations (6) 
15 3 OWDMA decoding derives estimates x(k|q) of x(k|q) for 

data block q from the receiver estimates z (i) of z(i) 
x(k|q) = FWT" 1 [ z(i) ] 

= Ei z(i) vj/*(i|q / k) 
4 FWT algorithm for OWDMA decoding in receiver 
20 4. a) FWT pre -sum 

X(iolq) = S A z (iolq+A) \|/(i 0 +AN) 
= S A £ (io+AN|q) V|/(i 0 +AN) 
4.b) FWT pre-sum FFT 

x(k|q) = Sio X(iolq) E(k / i 0 ) 
25 5 Computational complexity of fast algorithm 

Real multiply rate R M 

R M T = 2 log 2 (N) + 2 L 
Real add rate R a 

R A T = 3 log 2 (N) + 2 L 
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Step 4 is the new fast algorithm for the received OWDMA 
which partitions the baseband symbol detection x(k|q) in step 3 
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into a pre-sum calculation sub-step 4. a of Miolq) followed by a 
sub-step 4.b FFT of this pre-sum. 

Step 5 evaluates the real multiply complexity metric R M T 
and real add computational complexity metrics R A T in terms of 
5 multiplies/adds per digital sample for the fast algorithm in step 
3. The first term in these metrics is the complexity of the FFT' 1 
for a base 2 and the second term is the complexity of extending 
the multi-resolution complex Wavelet waveform over L of the 
symbol intervals. 

10 MS-CDMA parameters and codes are defined in equations (11) . 

Step 1 defines the scenario parameters. Step 2 partitions the 
user index u field into the sub-fields u 0 ,Ui of size N 0/ Ni for 
scales 0,1 respectively and which are the indices over the users 
within each channel and the indices over the channels within the 

15 MS-CDMA group and which uniquely represent u as u = u 0 +UiNi. 

Step 3 partitions the code chip index n field into the sub- 
fields n 0 ,ni of size N 0 ,Ni for scales 0,1 respectively and which 
are the indices over the chips within each channel and the 
indices over the channels of the MS-CDMA group and which uniquely 

20 represent n as n = n 0 +niNi. 

Step 4 defines the N c xN c MS-CDMA code matrix C whose 
elements are C(u,n) where u+1 is the row index and n+1 is the 
column index and where the +1 has been added to correspond to the 
row and column numbering starting with +1. MS-CDMA code vector 

25 

MS-CDMA parameters and codes (11) 
1 Scenario parameters 

M = number of communications channels 
of band B 

30 N 0 = number of CDMA chips per channel 

Ni = number of channels in MS-CDMA group 

N c = number of chips in MS-CDMA group = N 0 Nx 

1/T 0 = MS-CDMA chip/symbol rate 

= 1/NT for OFDMA, OWDMA in FIG. 1,2 
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x(u,q)= User u in group q 

= x(k|q) when N 0 =l, u=k 

= x(k) when N 0 =l , u=k, q=0 
User index u for a MS-CDMA group 

u = Index of CDMA users 
= 0,1,..., (Nc-1) 

= (u 0 ,Ui) field representation of u 
= u 0 + UiN 0 for fields u 0 , Ui 
where the index fields are 
Uo = index of users in a channel in 
the MS-CDMA group 
= 0,1,..., (No-1) 
Ui = index of channels in the MS-CDMA group 
= 0,1,...,N!-1 
Code chip index n for a MS -CDMA, group 
n = Index of CDMA code chips 
= 0,1,..., (Nc-1) 

= (n 0 ,ni) field representation of n 
= n 0 + niN 0 for fields n 0 , n x 
where the index fields are 

n 0 = index of code chips within a channel 
in the MS-CDMA group 
= 0,1,..., (N 0 -1) 
ni = index of channels in the MS-CDMA group 
= 0,1,..., ( Ni-1) 
MS-CDMA code matrix C for a MS-CDMA group 
C = N c x N c MS-CDMA code matrix 

= [C(u,n)] matrix of elements {C(u,n)} 
= [c(u)] matrix of lx N c code vectors 
C(u) = MS-CDMA code vector u which is row u+1 in 
C counting rows, columns starting with 1 
= [C(u,0), C(u,l), . . . ,C(u, N c -1)] 
Kronecker product example for C for a MS-CDMA group 

C = Ci®C 0 Kronecker product of Ci and C 0 
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= [ C (u,n)=Ci(ui,n 1 )Co(u 0 ,n 0 ) ] 



where 

® = Kronecker product of Ci and C 0 
Ci = Ni x Ni MS-CDMA scale 1 code matrix 
5 = [C!(ui,ni)] matrix of elements {Ci(ui,ni)]} 

over the user channels within the group 
C 0 = N 0 x N 0 MS-CDMA scale 0 code matrix 

= [Co(uo,n 0 )] matrix of elements {Ci(ui f rii)]} 
over the user chips within each channel 
10 in the group 

c(u) for user u is the row u+1 of C. The C is a complex 
orthogonal code matrix. 

Step 5 is an example of a MS-CDMA code matrix C which is 
15 factorable into a Kronecker product C=Ci®C 0 of C x and C 0 where 
is the Kronecker product and the matrix Ci is the Ni x Ni MS- 
CDMA scale 1 code matrix over the user channels within the MS- 
CDMA group and the matrix C 0 is the N 0 x N 0 MS-CDMA scale 0 code 
matrix for the user chips within each channel in the MS-CDMA 
20 group. 

A general definition for MS -CDMA, code matrix C which 
includes the definitions in references [3], [4] is allowed by this 
invention disclosure. The example Kronecker product construction 
of C is important because it is a property possessed by the real 

25 Walsh or Hadamard orthogonal matrices as well as the majority of 
orthogonal matrices used or proposed for CDMA and because of the 
relative ease in finding a fast algorithm for the transmitter and 
receiver coding and decoding and because of the ease in 
constructing multiple scales for MS-CDMA as demonstrated in 5 in 

30 equations (11) . 

MS -CDMA, representative application to OFDMA in FIG. 1 uses 
the candidate architecture which spreads the CDMA over the M 
contiguous 48 data channels or 52 contiguous data plus pilot 
channels and the representative application to OWDMA in FIG. 3 
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uses the candidate architecture which spreads the CDMA over the M 
contiguous OFDMA communications channels across B. 

MS-CDMA partitions the M channels into channel groups of 
size Ni and provides a scale over the chips within the channels 
5 and another scale over the channels within this group. Code chip 
length N c for each user in a channel group is equal to N c =N 0 Ni 
where N 0 is the number of chips in each channel assigned to scale 
"0" and Ni is the number of channels assigned to scale "1" 
within the group. Each user has a chip scale "0" and a channel 

10 scale "1" . Chip scale "0" spreads the data over the chips within 
each channel and channel scale "1" then spreads the channel chips 
uniformly over the Ni channels with the result that each user 
occupies each of the channels within the Ni channel group. 

There could be from 1 to M/Ni channel groups depending on 

15 the architecture and applications. The use of multiples groups 
M/Ni>l tends to be desirable since the storage requirements and 
computational complexity are reduced as the number of groups are 
increased and the spreading advantages within each group tend to 
saturate as the number Nx>16 when the channels within each goup 

20 are spread across the fullband M channels. 

For N 0 =l there is no CDMA within each channel and the MS- 
CDMA then spreads the signals over each channel within a group 
for both OFDMA and OWDMA to function as a means to spread each 
channel over the fullband M channels and which may be a desirable 

25 architecture when the channels are sufficiently narrow to produce 
a sufficiently long pulse to counter multipath. 

MS -CDMA OFDMA transmitter equations are defined in 
equations (12) for MS-CDMA. Step 1 gives the parameters and 
definitions. Step 2 defines the encoding equations for chip n 0 

30 for data block q. Sub-step 2. a uses the fast code transform FCT 
developed in references [3] , [4] to generate the encoded vector. 
Sub-step 2.b uses the inverse FFT" 1 to construct the transmitter 
complex baseband signal z(i 0 |n 0 +qN 0 ) for chip n 0 for data block q 
and these signals are combined in sub-step 2.c to generate the 

35 transmitter signal z(i) for all n 0 ,q. 
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Step 3 evaluates the real multiply complexity metric R M T 
and real add computational complexity metrics R A T in terms of 
multiplies/adds per digital sample for the fast algorithm in step 
3. The first term in these metrics is the complexity of the FFT" 1 

5 

MS -CDMA OFDMA transmitter equations (12) 

1 Parameters and definitions are defined in 1 in 
equations ( 2), (9), (11) together with 

MS-CDMA group index field {g} 
10 g = 0,1,. - - , (M/Ni-1) where M/N x = number 

of MS-CDMA groups over the M data channels 
Frequency harmonic k 

k = k(nilg,q) 
MS-CDMA chip index fields 
15 n = n 0 + niN 0 

= 0,1,2, . . ,N C -1 
n 0 = 0,1,2,. . ,N 0 -1 scale 0 chips in channel 
ni = 0,1,2,. . ,Ni-l scale 1 chips over channels 
N c = N 0 Ni number of chips in a MS-CDMA group 
20 MS-CDMA channel code index ni 

ni = ni(k|g,q) function of k for a given g,q 
User symbol x 
x = x(u|g,q) 

= A(u|g,q) exp j(p(u|g,q) 
25 MS-CDMS code orthogonality 

CC* = N c I 

PN spreading code P(n) for chip n in group g 

P(n|g,q) = exp ( jcp(n|g,q) 
Assume the band and pulse filtering can be 
30 neglected on both Tx and Rx which is 

equivalent to the assumption p(i)=l for all i 

2 MS-CDMA OFDMA. encoding f or complex baseband signal z (i) 
2. a) FCT MS-CDMA encoding for data block q 

Z(n|g,q) = FCT [ x(u|g,q) ] 
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= N c * x £ u x(u|g,q)C(u,n) P(n|g,q) 
where FCT is the fast code transform 
2.b) FFT" 1 OFDMA encoding for chip n 0 of data block q 
z(i 0 |n 0 +qN 0 ) = FFT' 1 [ Z(n|g,q) ] 
5 m N' 1 Z k Z(n|g,q)E*(k,i 0 ) 

where n = n 0 + niN 0 
ni = ni(k|g f q) 
Z (n| g,q) =Z (n 0 +n! ( k | g, q) N 0 ) 
2.c) Complex baseband transmitter signal z (i) 

10 z(i) = Z( 0 , z(i 0 |n 0 +qN 0 ) for (o) = (n 0 +qN 0 ) N 

3 Computational complexity 

Real multiply rate R M 

R M T = 2 log 2 (N) 
Real add rate R a 
15 R A T = 3 log 2 (N) + 2 log 2 (N c ) + 2 

for a base 2 and the second term is the complexity of the FCT 
assuming the FCT does not require any multiplies. 

MS-CDMA OWDMA transmitter equations are defined in 

20 equations (13) for MS-CDMA. Step 1 lists the parameters and 
definitions. Step 2 defines the encoding equations for chip n 0 
for data block q. The FCT on the symbols in sub-step 2. a yields 
the encoded data block Z(n|g,q) and the FWT on this output in 
sub-step 2.b yields the transmitter complex baseband signal 

25 z(i 0 |n 0 +qN 0 ) for chip n 0 for data block q and these signals are 
combined in sub-step 2.c to generate the transmitter signal z(i) 
for all n 0 ,q. 

MS -CDMA. OWDMA transmitter equations (13) 
30 1 Parameters and definitions are defined in 1 in 

equations (9), (12) together with the range of q 
q = 0 for the data block being addressed 
= +/-1 contiguous data blocks 
Assume the h(f) is flat over the passband for 
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both Tx and Rx and can be neglected 
2 MS -CDMA OWDMA encoding for chip n 0 for data block q 
2. a) FCT MS-CDMA encoding for data block q 

Z(n|g,q) = FCT [ x(u|g,q) ] 
5 = Nc" 1 S u x(u|g,q)C(u,n)P(n|g) 

2.b) FWT OWDMA encoding for chip n 0 of data block q 

z(i 0 |n 0 +qN 0 ) = E 8q Z A FWT [ Z (n+A| g, q+8q) ] 
= N" 1 Z k E 5q 2 A Z (n+A|g,q+5q)i|/(i |n 0 +A,k) 
= N^ZkSaqSA Z (n+A|g,q+8q)i|/(i 0 +AN)E* (k,i 0 ) 
10 where the overlap conditions on 5q are 

8q = A for N 0 =l 

= sign (A+n 0 ) L I A+n 0 1 /N 0 J otherwise 
sign(o) = numerical sign of (o) 
L(o)J = rounded down integer value of (o) 
15 where the overlap index A = 0, +/-1,. 

. . .,+/-(L/2-l) 
2.c) Complex baseband transmitter signal z(i) 

z(i) = I (0) z(i 0 |n 0 +qN 0 ) for (o) = (n 0 +qN 0 ) N 

3 FWT algorithm for OWDMA encoding in the transmitter 
20 3. a) FWT pre-calculation FFT" 1 

Mioln 0 +A,q) 

= N" 1 E k Z(n+A|g,q)E*(k,i 0 )] 

= FFT" 1 [IT 1 I k Z(n+A|g / q)E*(k / i 0 )]] 
3.b) FWT post-sum 

25 z(i 0 |n 0 ) = SqS A v(i-(n 0 +A)N)X(ioln 0 +A # q) 

4 Computational complexity for fast algorithms in 2. a, 3 

Real multiply rate R M 

R M T = 2 log 2 (N) + 2 L 
Real add rate R a 
30 R A T = 3 log 2 (N) + 2 L + 2 log 2 (N c ) + 2 
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Step 3 is the new fast FWT algorithm in this invention 
disclosure for the transmitted OWDMA which consists of the pre- 
calculation FFT" 1 in sub-step 3. a followed a post-sum in sub-step 
3.b of the product from sub-step 3. a with the corresponding 
5 Wavelet overlaps over the q data block interval. 

Step 4 evaluates the real multiply complexity metric R M T 
and real add computational complexity metrics R A T in terms of 
multiplies/adds per digital sample for the fast algorithms FCT 
and FWT in steps 2,3. The first term in these metrics is the 

10 complexity of the FFT" 1 for a base 2, the second term 2L is the 
complexity of extending the multi-resolution complex Wavelet 
waveform over L of the MT symbol intervals, and the remaining 
terms are the complexity of the FCT assuming the FCT does not 
require any multiplies. 

15 MS-CDMA OFDMA/OWDMA encoding for the transmitter in FIG. 6 

is a representative implementation of the MS-CDMA OFDMA encoding 
algorithms in equations (12) and the MS-CDMA OWDMA encoding 
algorithms in equations (13) . Signal processing starts with the 
input stream of data encoded symbols x(u|g,q) 40 from the 

20 transmitter symbol encoder 52 in FIG. 7 and defined in 1 in 
equations (12) for both OFDMA and OWDMA. 

FIG. 6 MS-CDMA encoding 41 implements the fast cpde 
trams form FCT encoding defined in sub-step 2. a in equations 
(12), (13) for the MS-CDMA encoding and PN cover or spreading 

25 encoding to generate the N c outputs x (u I g, q) C (u, n) P (n | g) for each 
MS-CDMA group g and data block q. For each group these outputs 
are summed 42 over u to generate the encoded vector Z(n|g,q) in 
43. 

OFDMA processing 44 performs an FFT" 1 on the received set of 
30 vectors Z(n|g,q) and a summation to implement sub-steps 2.b,2.c 
in equations (12) and the output is band filtered 46 to generate 
the MS-CDMA OWDMA encoded complex baseband signal z(i) in 47. 

OWDMA processing 45 performs an FWT on the received set of 
vectors Z(n|g,q) and a summation to implement sub-steps 2.b,2.c 
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in equations (13) and the output is band filtered 46 to generate 
the MS-CDMA OWDMA encoded complex baseband signal z(i) in 47. 

Outputs z(iolno) 47 from the MS-CDMA OFDMA and MS-CDMA OWDMA 
are digital-to-analog DAC converted 48 and handed off to the 
5 analog front end 49 as the complex baseband analog signal z(t) in 
49. 

MS -CDMA. OFDMA/ OWDMA transmitter block diagram in FIG. 7 is 

a representative transmitter implementation of the MS-CDMA OFDMA 
and MS-CDMA OWDMA encoding in the shaded box in FIG. 6. The 

10 transmitter block diagram in FIG. 7 includes the FIG. 6 MS-CDMA 
OFDMA/OWDMA encoding in an abbreviated format in the shaded box 
54 in FIG. 7. FIG. 7 signal processing starts with the stream of 
user input data words. Frame processor 51 accepts these data 
words and performs the encoding and frame formatting, and passes 

15 the outputs to the symbol encoder 52 which encodes the frame 
symbols into amplitude and phase coded symbols x(u|g,q) 53 which 
are the input to the MS-CDMA encoding 55 and which is 41,42 in 
FIG. 6. 

MS-CDMA FCT encoding outputs Z(n|g,q) 56 are handed over to 
20 the OFDMA and OWDMA processing 57 which performs an inverse FFT" 1 
followed by a band filtering for OFDMA which is 44,46 in FIG. 4 
and performs an inverse FWT" 1 followed by a band filtering which 
is 45,46 in FIG. 6. This complex baseband signal z(i) in 47 in 
FIG. 6 is digital-to-analog DAC converted 59 and the the output 
25 complex baseband analog signal z(t) 60 is handed off to the 
analog front end 61. 

The z(t) is single sideband upconverted, amplified, and 
transmitted (Tx) by the analog front end 61 as the real 
waveform v(t) 62 at the carrier frequency f 0 whose amplitude 
30 is the real part of the complex envelope of the baseband waveform 
z(t) and the phase angle $ accounts for the phase change from the 
baseband signal to the transmitted signal. Output waveform 62 
from the analog front end is the Tx waveform from the Tx antenna. 
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It should be obvious to anyone skilled in the 
communications art that this example implementation in FIG. 
6,7 clearly defines the fundamental MS-CDMA OF DMA and MS-CDMA 
OWDMA signal processing relevant to this invention disclosure 
5 and it is obvious that this example is representative of the 
other possible signal processing approaches. 

MS -CDMA OFDMA receiver equations are defined in 
equations (14) , Step 1 lists the parameters and definitions and 
the assumption that the band and pulse filtering can be 
10 neglected. Step 2 defines decoding of the received chip signal 
to derive the estimate x(u|g,q) of the transmitted symbol 
x(u|g,q). Sub-step 2. a derives the estimate for the encoded 
symbols using the FFT and sub-step 2.b uses the FCT" 1 on this 
estimate to derive the transmitted symbol estimate. 

15 

MS -CDMA OFDMA receiver equations (14) 

1 Parameters and definitions are defined in 1 in 
equations (2), (9), (11) together with 

z (i) = receiver band filter estimate of 
20 the transmitted complex baseband signal z(i) 

Assume the band and pulse filtering can be 
neglected on both Tx and Rx which is 
equivalent to the assumption p(i)=l for all i 

2 MS-CDMA OFDMA decoding to derive data symbol x (u|g,q) 
25 2. a) FFT of each received chip vector z (i 0 |n 0 +qN 0 ) 

Z(n|g,q) = estimate of transmitted Z(n|g.q) 
- FFT [ z (ioln 0 +qN 0 ) ] 
= Z i0 z (i 0 |n 0 +qNo) E(k,i 0 ) 
where z(i 0 |n 0 +qN 0 ) = z (i=i 0 + (n 0 +qN 0 ) N ) 
30 2.b) FCT" 1 of each group g encoded vector Z(n|g,q) 

x (u|g,q) = FCT" 1 [ Z (n|g,q-0) ] 

Z n Z (n|g,q)P*(n|g)C*(u,n) 
3 Computational complexity is the same as calculated 
in 3 in equations (12) for the transmitter 
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MS-CDMA OWDMA receiver equations are defined in equations 
(15) . Step 1 lists the parameters and definitions and the 
assumption that the band and pulse filtering can be neglected. 
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MS-CDMA OWDMA receiver equations (15) 

1 Parameters and definitions are defined in 1 in 
equations (2) , (9) , (11) , (14) together with 

10 z (i) = receiver band filter estimate of 

the transmitted complex baseband signal z(i) 
Assume the h(f) is flat over the passband for 
both Tx and Rx and can be neglected 

2 MS -CDMA OWDMA decoding to derive data symbol x (u|g,q) 
15 2. a) FWTT 1 of each received signal z (i) 

Z(n|g,q) = estimate of transmitted Z(n|g.q) 
= FWT _1 [ z (i) ] 
= S ± z (i)\|/*(i|q f k) 

= Z i0 2 A z (i 0 |n 0 +(q+A)N)v)/(io+AN)E(k / io) 
20 where z(i 0 |n 0 +qN 0 ) = z (i=i 0 + (n 0 +qN 0 ) N ) 

k -> g which means k specifies g 
n 0 ,k -> n=n 0 +ni(k|g,q) 
2.b) FCT" 1 of each group g encoded vector Z(n|g,q) 

x (u|g,q) = FCT" 1 [ Z (n|g f q) ] 
25 Z n Z (n|g,q)P*(n|g)C*(u,n) 

3 FWT" 1 algorithm 

2. a) FWT" 1 pre-s\am 

^(iolg.q) = S A z (i 0 |n 0 +(q+A)N)v)/(i-(no+A)N) 
where n 0 +A g using the boundary conditions 
30 in 2.b in equations (13) 

2.b) FWT" 1 pre -sum FFT 

Z (n|g # q) = S i0 FFT [ ^(i 0 lg,q) ] 
= Z i0 X(i 0 lg/q)E(k,i 0 ) 
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3 Computational complexity is the same as calculated 
in 3 in equations (13) for the transmitter 



Step 2 defines decoding of the received chip signal to derive 
5 the estimate x(u|g,q) of the transmitted symbol x(u|g,q). Sub- 
step 2. a derives the estimate for the encoded symbols using the 
FWT" 1 and sub-step 2.b uses the FCT" 1 on this estimate to derive 
the transmitted symbol estimate. Step 3 is the new fast FWT 
algorithm in this invention disclosure for the received MS-CDMA 

10 OWDMA which consists of the pre-sum in sub-step 3. a followed an 
FFT on this pre-sum. 

MS -CDMA. OFDMA/OWDMA decoding for the receiver in FIG. 8 is 
a representative implementation of the MS-CDMA OFDMA decoding 
algorithms in equations (14) and the MS-CDMA OWDMA decoding 

15 algorithms in equations (15) . Signal processing starts with the 
input intermediate frequency IF signal after being single 
sideband downconverted and synchronized in frequency 63 from 
the receiver front end 72 in FIG. 9. This input signal is band 
filtered 64 and handed off to the analog-to-digital converter ADC 

20 or A/D 65 whose digital output is the received estimate z(i) 66 
of the transmitted complex baseband signal z(i) 60 in FIG. 7. 
For OFDMA 67 this received signal is processed by an FFT to 

derive estimates Z(n|g,q) 68 of the MS-CDMA encoded signal 
Z(n|g,q) 56 in FIG. 7. For OWDMA 67 this received signal is 
25 processed by an inverse multi-resolution complex Wavelet 

transform FWT* 1 to derive estimates Z(n|g,q) 68 of the MS-CDMA 
encoded signal Z(n|g,q) 56 in FIG. 7 implementing the fast 
algorithm defined in 2 in equations (15) . Recovered estimates 

Z(n|g,q) 68 are processed by the inverse fast code transform 
30 FCT* 1 69 to derive estimates x(u|g,q) 70 of the transmitted data 
symbols x(u|g,q) 53 in FIG. 7 for hand over to the symbol 
decoder. 

MS-CDMA OFDMA/OWDMA receiver block diagram in FIG. 9 is a 

representative receiver implementation of the MS-CDMA OFDMA and 
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MS-CDMA OWDMA decoding in the shaded box in FIG. 8. The receiver 
block diagram in FIG. 9 includes the FIG. 8 MS-CDMA OFDMA/ OWDMA 
decoding in an abbreviated format in the shaded box 74 in FIG. 
9. FIG. 9 signal processing starts with the received Rx waveform 
5 71 from the transmitter 62 in FIG. 7. Received signal v(t) 71 
is an estimate of the transmitted signal v(t) 62 in FIG. 7 
received with errors in time At, frequency Af, phase AG, and 
with an estimate z (i) of the transmitted complex baseband signal 
z(t) 60 in FIG. 7. This received signal v(t) is amplified and 

10 downconverted by the analog front end 72 and single side band 
SSB downconverted and synchronized 73 and handed over to the 
digital-to-analog conversion DAC processing 75 for band filtering 
by h(f) and digitization to generate the baseband signal z(i) 76 
which is the received estimate of the transmitted signal z(i) 60 

15 in FIG. 7. Timing synchronization could be implemented in the 
DAC. 

Outputs z (i) are processed by the MS-CDMA OFDMA/ OWDMA 
decoding to derive estimates x(u|g,q) 79 of the transmitted 
symbols x(u}g,q) 53 in FIG. 7 and part of the information is 
20 handed off to the synchronization (sync) processor. For the 
inverse MS-CDMA" 1 OFDMA" 1 the processing 79 consists of the FFT 

to recover estimates Z(n|g,q) of Z(n|g,q) followed by an inverse 
FCT" 1 to recover x(u|g,q). For the inverse MS-CDMA" 1 OWDMA" 1 the 
processing 79 consists of the inverse FWT" 1 to recover estimates 

25 Z(n|g,q) of Z(n|g,q) followed by an inverse FCT' 1 to recover 
x(u|g,q) and the FWT" 1 implementes the fast transform in 2 in 
equations (15). Outputs x(u|g,q) are processed by the symbol 
decoder 80 and the frame processor 81 for handoff as the 
received Rx data 82. 

30 It should be obvious to anyone skilled in the 

communications art that this example implementation in FIG. 
8,9 clearly defines the fundamental MS-CDMA OFDMA and MS-CDMA 
OWDMA signal processing relevant to this invention disclosure 
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and it is obvious that this example is representative of the 
other possible signal processing approaches. 

Variable power control across the frequency band can be 
implemented by assigning each group g of transmit Tx signals 
5 their own power level P(g). Each MS-CDMA group g occupies a 
subset of the channels over the frequency band B consisting of 
Ni channels which means that the users within group g are 
transmitted with the same Tx power. On receive each group g of 
channels is processed separately so there is no cross-talk 

10 between the users in the different groups. The OWDMA was designed 
to support large dynamic range imbalances between channels which 
could be present with power level control. MS-CDMA OFDMA/OWDMA 
variable power control over the frequency subbands corresponding 
to the MS-CDMA groups supports the potential for diversity 

15 improvements by allocation of the available power to emphasize 
the *best' set of available subbands, which subbands are not 
required to be contiguous as well as supporting the simultaneous 
support of multiple users with differing power requirements due 
to range, multi-path, and path attenuation effects. 

20 A second configuration for variable power control is 

described in reference [2] and increases the flexibility of power 
control to all of the individual channels. 

Preferred embodiments in the previous description is 
provided to enable any person skilled in the art to make or use 

25 the present invention. The various modifications to these 
embodiments will be readily apparent to those skilled in the art, 
and the generic principles defined herein may be applied to other 
embodiments without the use of the inventive faculty. Thus, the 
present invention is not intended to be limited to the 

30 embodiments shown herein but is not to be accorded the wider 
scope consistent with the principles and novel features disclosed 
herein. 
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